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ABSTRACT 

Waves reflected by the inner edge of a protoplanetary disk are shown to 
significantly modify Type I migration, even allowing the trapping of planets near 
the inner disk edge for small planets in a range of disk parameters. This may 
inform the distribution planets close to their central stars, as observed recently 
by the Kepler mission. 

Subject headings: hydrodynamics — planet-disk interactions — planets and satel- 
lites: dynamical evolution and stability — protoplanetary disks — waves 
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1. Introduction 

The location and configuration of planets within their planetary systems depend sen- 
sitively on the effe cts of many physical pro cesses including planet-plan et scattering (e.g. 
Rasio fc Fordlll996l). Kozai oscillation s (e.g. iFabryckv &: Tremaind 120071 ). and disk migra- 
tion (e.g. iGoldreich fc TremaindllQSOl . hereafter GT80)lj. Of these, disk migrati on is likely to 



have been dominant in systems where planets are observed to be co-planar (e.g. iSteffen et al. 
2010l ). The density waves launched inwards from the inner Lindblad resonances carry away 
negative angular momentum driving the planet outwards (positive torque). The waves 
launched outwards from the outer Lindblad resonances carry away positive angular mo- 
mentum driving it inwards (negative torque). For planets too small to clear a gap in their 
protopla netary disk Type I migration occurs due to the asymmetry between the Lindblad 
torques (IWardI Il986l . Il997l . hereafter W86 and W97) with the outer torque being typically 
stronger, driving a net migration inwards. 

Recent observations by the Kepler spacecraft have shed new light on to t he distribution 
of extrasolar planets within their planetary systems. [Howard et al.l (1201 ll ) examined the 



^For review and comparison of migration mechanisms see iMorton fc JohnsonI (J2011I ) and lWard fc Hahn 



(J2000l ). and references therein. 
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radial distribution of planets closer than 0.25 AU of solar type stars observed in the Kepler 
field. They found that smaller planets {Rp < 2 — 4/?^) have an exponential cutoff in the 
radial distribution at an orbital period ~ 7 days, while larger planets have exponential cutoffs 
closer to the star at ~ 2 days. 

The protoplanetary disks are likely be truncated before direct contact wit h the host 



stars by a strong stellar magnetic field near to the corotation radius (see e.g. IShu et al. 
19941 ). Rotational studies of classical T Tauri s tars give periods ranging from ~ 3 — 7 days 



for stars with M* > O.25M0 (JLamm et al.ll2005l ). The truncation radius is typically taken to 
be ~ 5 stellar radii, or abou t ~ 4 day periods for typical 2 solar radius pre-main-sequence 
stars (JGuUbring et al.lll998[ ). The inner edge of the accretion disk plays a role in halting 
migration very close to the ho st star; disk mig ration is naturally halted once a planet passes 
through the inner disk edge (JLin et al.l Il996l ) and the Lindblad resonances can no longer 
interact with the disk material. 

Close to the inner edge of the disk the effect of the corotation resonance also becomes 
important and can exceed the differential Li ndblad torque, particularly at locations with 
strong gradients in the surface density profile ( iGoldreich fc Tremaindl 19791 ). However, unless 
the viscosity is strong enough to transfer angular momentum away from the co-orbital region 
quickly enough, the corotation resonance saturates and the cor otation torque goes to zero 
jMassetll200l[IOgilvie fc Lubowll2003l : ISaruteau fc Massed [2008[ ). 



In this letter we will outline a halting mechanism for Type I migration of planets due 
to density waves resonant between the inner Lindblad resonance and the disk inner edge. 
This allows for the Type I migration of smaller planets to be halted within the disk, at a 
few times the inner disk radius. 



2. Torque due to Lindblad Resonances 

2.1. Basic Equations 

To evaluate the effect of a perturbing body on a gas disk we limit ourselves to the 
examination of the 2-dimensional disk and proceed with the vertically integrated disk vari- 
ables S = f p dz and P = J p dz the surface density and vertically integrated pressure 
respectively. We assume that in the unperturbed state the disk velocity profile is Keplerian, 
Uq = rQ(f) and that perturbations are barotropic such that P = -P(S). 



We begin by first considering the inviscid limit of the system. The linear perturbation 
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equations for continuity and momentum conservation are 



dt 



+ V ■ (S5u + Uo5S) = , (1) 



^ + (u, ■ V)Su + {6u ■ V)u„) = -^V6P + gvP - V$, (2) 

where 6 signifies the hnear perturbation of a quantity, $ is the external potential due to the 
perturbing planet with mass Mp located at radius Vp away from the central star. 

For a perturber with vanishing eccentricity we can break dow n this perturbing potential 



into azimuthal Fourier components (JGoldreich fc Tremainelll979l . hereafter GT79) 



$ = y, ^m{r) exp(^m0 — imVtpt) , (3) 

m=0 

where Vtp is the orbital frequency of the perturber, and 

rp ^2 ^-''> [by,{r/rp) ^^^ 
Si J is the Kronecker delta function and 6^2 (^) is the Laplace coefficient: 

cos mOdO 
ir Jo (l-2a;cos^ + a;2)V2 



^rn = -^(2 - 6rr,o) [bTMr/Tp] - ^^n.,! ) , (4) 



bTd^) - Z / .. ..,.„.., .,.^W. • (5) 



For any given Fourier component $„, the disk response will have perturbations of the 
form 6 oc exp{im(j) — irriQpt) giving 

/ 1 f/ ITYl 

— iCj—5h H — -—-{rT,5ur) H T^dus = 0, (6) 

cj. r or r 

d 

— iuj5ur — 2Qs5us = — 7— (5/i + $„), (7) 

or 
2 
K im 

— iu6u^ + —5ur = {5h + ^rn), (8) 

where k^ = ^^(r^f2) is the square of the radial epicyclic frequency, u = m{Qp — Q) is 
the Doppler shifted perturbation frequency, 6h = 6P/T, is the enthalpy perturbation and 
Cg = dP/dT, is the square of the sound speed. 

Eliminating 6ur and ^m^ we find 

C{5h + ^^) + ^$„, = 0, (9) 
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where the operator C is defined as 



rf2 /^ ^Y.\ d 
dr"^ \dr D J dr 



2mn f d , VtT. 

— ^ ~r ^^ — 
ruo \dr D 



m 



D 

72 



,. ,. (10) 

and D = k^ — oJ ^- This is our basic working equation (equivalent to equation (29) in 
Zhang fc Laill2006l . hereafter ZL06, with n = 0). 



2.2. The Differential Lindblad Torque 

Following §4 in ZL06 we can expand Eq. i^ in the vicinity of a Lindblad resonance r^ 
where D{rL) = 0, in the cold disk approximation. Changing variables to x = (r — tlj/tl 
and keeping the important singular terms this gives 



dx^ X dx / X 



where 



and 



^. 



d 2mn\ 
r- z— 1 $r 

dr u 



/3 = 



1 dD 



cldhir 



:in 



(12) 



(13) 



rL 



To solve the above equation we define a new variable t] such that 5h = drj/dx. Integrating 
and rearranging, Eq. (ITT]) becomes 



dx'^ 



rj — (5xr] = — \&^ 



(14) 



where we can absorb any integration constant into the definition of t]. 



The homogeneous solutions for Eq. (1141) are the Airy functions (Abramowitz & Stegun 
1964) 

r/i = Ai{(3'/'x), (15) 

r/2 = Bi(/3i/\), (16) 

while the general solution can be found using the method of variation of parameters 



(17) 



^ = TT7 I ^72 / Vl'^mdx -Tji I TJi"^ mdx 

+My^ + Ny2 . 
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where W = rii{dri2/dx) — 772(^771 /dx) = fi^l"^ jix is the Wronskian, and M and A^ are constants 
determined by the boundary conditions. 

If the function \l/m does not vary strongly over the first Airy half-wavelength near the 
resonance then the solution can be approximated by 



^ - ^2/3 



Ai(Orfe + iV)Bi(0 


Bi(Oc?e + M)Ai(0 





(18) 



where i = P^l^x, M = M/32/3/(7r^^(ri)) and A^ = iV/32/3/(7r^^(ri)) 



The asymptotic expansions of the Airy functions are well known. For .^ < the behavior 
of the Airy functions is wavelike, while for ^ > the Ai(^) behaves as a decaying exponential, 
while Bi(^) grows exponentially. Since Bi(^) grows exponentially and the physical solution 
must be bounded we require that as ^ — ;■ 00 the term multiplying Bi(^) must vanish. Hence 
we find that 

N = - P kmdi = -\. (19) 

To determine M we must consider the wave region in detail. The asymptotic expansions 
for Ai(^) and Bi(^) for ^ <^ — 1 are (taken from ZL06) 

Ai(0 ~ 7r-i/2(-0-^/'sinX(0, (20) 

Ai(0 ~ -^ + 7r-^/2(-0"'/'cosX(0, (21) 

Bi(0 ~ n~''\-ir'"co^X{0. (22) 

^Bi(0 ~ -7r-V2(-0"'/'sinX(0, (23) 



X{^) = l{-e" + \- (24) 



where 



This gives the asymptotic solution 



^h-'^^{-if'% 1 -.M)e^^«) 



+ (l + ^M)e-*-^^«^]. (25) 

As in ZL06 we have for the outer Lindblad resonance (roLn) an outgoing wave if the 
solution has the form exp(-|-«X(^)). This outgoing wave solution is achieved by setting 
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M = i, and yields 



Sh{r > tolr) 



/31/3 



-0 



l/4giX{0 



(26) 



For the inner Lindblad resonance (riLn) an ingoing wave (propagating away from the 
resonance) has the form exp(— zX(,^)) (since the group velocity and phase velocity have 
opposite sign inside the corotation resonance), and an outgoing wave (propagating towards 
the resonance) has the form exp(+iX(^)). With no reflecting inner boundary in the disk the 
ingoing wave correctly satisfies the boundary condition (no waves are reflected back) and 
the correct solution is obtained by setting M = —i such that 



6h{r < tilr) 



' /3V3 



-0 



l/4g-iX(0 



(27) 



The torque can be calculated by determinin g the angular momentum flux (see GT79; 
Lvnden-Bell fc Kalnaisl[l972[ IXanaka et al.ll2002[ ) 



F(r) 



vrmrE / d6h* 
-Im oh- 



D 



dr 



(28) 



where * denotes the complex conjugate. The torque on the disk due to each resonance is 
then given by the flux evaluated in the wave region away from the resonance. This gives 



Tm(roLR) = F{r > tolr) 



mn^E'^l 



for the outer Linbdlad resonance and 



Tm{riLK) = F{r < hlr) 



rdD/dr 



mvr^Svl/^ 



(29) 



^OLR 



HLR 



rdD/dr 
for the inner Lindblad resonance, which agrees with W86 and GT79. 



(30) 



As discussed in W97 and lArtymowica (Il993l . [l99J) the torque cutoff/pressure buffer 
effect can be included if \E'm is multiplied by a factor of [1 + 4:(mc/ {r k,)Y]~^''^ , the D in ( l29l) 
and (l30l) is replaced by D' = k"^ — u"^ + m'^c^/r'^ (to account for the excess pressure due to 
azimuthal winding), and the torque expressions are instead evaluated at the modified inner 
and outer Lindblad resonances where D' = 0. 

Evaluating the inner and outer Lindblad torques we see in Figure [T] that for waves 
traveling away from resonances the outer torque will be stronger than the inner torque and 
the net migration must therefore be inwards in Type I migration. 
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Fig. 1. — The normal inner and outer Lindblad torques for various values of the disk thickness 
H/r evaluated including the effect of the pressure buffer (see W97), and normalized by the 
value To = {Mp/M^:)'^'Er'^Q'^{H/r)~^ evaluated at the planet location, where M* is the mass 
of the host star. The torques are cutoff for m > Amcutofr oc {H/r)~^. 



2.3. Torque Due to a Trapped Wave 



We can generalize the above results (see iTsang fc Laill2008l . for a more formal discussion 



for the homogeneous solution) for waves traveling far from the resonance by defining 

k ~ VL'/c2 (31) 

and substituting 

X(0 = ± r kdr + J, (-0^/^ = TrLk//3'/-' , (32) 

where the signs are for the outer and inner Lindblad resonances respectively. 

If there exists a boundary at some position rin < Hlr with a complex reflection coefficient 
TZ then 

5/i(rin) oc exp ii k dr \ + TZexp I — W k dr \ . (33) 

To satisfy Eq. f p3|) and (|25|) we must then have 

1 -iM 



1 + iM 



Tie'''' , (34) 



where 



This gives us 



M 



e= kdr + 71/4:. (35) 



g[J0+(ln7e)/2] _ g-[ie+(ln7^)/2] 



g[ie+(ln7^)/2] _^ g-[ie+(ln7^)/2] 

-tan[e-^(ln7^)/2]. (36) 



-mn^iM) ( f^"\ ) . (37) 



Applying Eq (ES]) to Eq. ([28]) we find 

Tm(riLR) = F{r < riLR 

rdD / dr J ^, 

When 7?. = we have M = —i recovering the no-refiection case. However for non-zero 
7^ we have 

M = - tan[e - Arg(7^"l/2) + i In |7^^^/2] . (38) 
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When O — Arg(7^ ' ) = Tr(n + 1/2) and the waves launched at the Lindblad resonance 
becomes trapped then we have M = ~ tanhffiii[R|-V2')] which gives the torque 

_ ( mn^J:K \ ( 1 \ 

V rdD/dr J^^^^ Vtanh[ln|7^|-V2]y " ^""^^ 

As \TZ\ —7- 1 (e.g. for the free boundary condition 6h = we have R = e^'^) the torque due 
to the resonant wave approaches T^ — ?■ oo. However, we must also take into account the 
damping of the propagating wave. 



2.4. Wave Damping 

If the damping is too strong, such that the damping length is comparable to the first 
Airy half- wavelength, then the torque cannot be easily related to the flux as in Eq. fl28|) . 
However, for damping length much greater than the driving length we may continue to use 
the flux evaluated near the Lindblad resonance as a measure of the net torque exerted on 
the disk. 

Introducing a small damping term into the system will contribute an imaginary compo- 
nent to the wavenumber, k = kr + iki, such that 

e = e^ + «6i = / kr + ikidr , (40) 



which gives the result 

— — -^ t tan[e, + ^e, + ^ In n^'^'] . (41) 

rdD/dr J^^^^ 

When Qr + Arg(7?.~^/^) = 7r(n -|- 1/2) and the density wave is trapped between the inner 
boundary and the inner Lindblad resonance we have the torque 

_ ( mn^^^l \ ( 1 \ 

yrdD/dr )^^^^\tanh[ln \n\~y^ + Qi] J ' ^ ' 

We define the resonant gain F^ as 

T 1 

^ (TmU^o tanh[ln|7^|-l/2 + e^] ' ^ ^ 

the gain in the inner Lindblad torque due to the wave being resonant in the disk. 
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Here, for simplicity, we take the viscosity to be independent of the perturbed quantities 
so that viscous damping takes the form of given by W97 

^^JSL^ll^I^, (44) 



where A ~ a ^m ^'^(H/r), H = Cs/Vt is the disk scale height, and a is the lShakura fc Sunyaev 



(119731 ) viscosity parameter. 



For strong wave damping 6j is large and Tm — ^ 1, the normal Lindblad torque. However, 
for weak damping 6j — > with perfect reflection \TZ\ — t- 1 we have F^ — ;• C)0, allowing for 
very large inner Lindblad torque.. 



3. Halting Type I Migration: Linear Case 

For the cases where the corotation torque is negligible. Type I migration will be halted 
if the total torque due to the inner Lindblad resonances exceed the torques due to the outer 
Lindblad resonances. The resonant gain as a function of Vp/ri^ are shown for the first few m 
in Figure [21 

There exists a torque cutoff (GT80, W97) such that only a finite range of ra contribute 
strongly to the total torque (see Figured]). As the planet migrates through the disk it sweeps 
through various resonances for the global modes. If the excess torque generated by a ??7,-mode 
that is in resonance exceeds the torque deficit due to all the other m within the torque cutoff 
then the migration can be stopped by the trapped wave resonance. 

We can estimate the maximum viscosity parameter a such that the migration can be 
halted by this effect by setting the maximal gain F^ roughly equal to the width of the torque 
cutoff 

Am^utoff = f (-^) , (45) 

where / is a factor of few that absorbs the variation of the torque for azimuthal modes 
within the torque cutoff. This is the worst case assumption, applicable since away from the 
strong resonant peak where the gain is F^es = l/tanh(6i + In |7^|^/^), we have a reduced 
gain, ~ l/Fj-es <^ 1- Thus the azimuthal mode in resonance must make up for, at most, the 
rest of the modes within the cutoff. We then have 

If we take the boundary condition to be 6h{rin) = then we have \TZ\ = 1, and 



Qi <, tanh 



f\r 



1 (V 1 ■ 
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Fig. 2. — Resonant gain Tm for the first three azimuthal modes for hnear viscous damping 
with H/r = 0.05 and a = 10~^. Far from the disk edge F^ — ?■ 1- The gain increases strongly 
for waves that are resonant in the disk. 
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where the last approximation holds for small H /r. The above condition says that Type 
I migration can be halted, if the number of e-foldings of wave amplitude lost across the 
resonant cavity is smaller than l/A?Ticutoff- 

This yields, for viscous dissipation. 

If we assume that rjn is located at the minimum distance for a resonance, the location of the 
first node, then we have the wavelength condition 

dr + 7r/4 ~ 7r/2 . (49) 

Near the Lindblad resonance we have D c:^ {r — rii^^{dD / dr)r^-^^ giving 

-1/2 

(50) 



/ \ 3/2 o 



6m 



Tp J r.p 



V H 
Thus if the viscosity parameter a is smaller than 

type I migration can be halted by trapped modes in the disk. 

The net Lindblad torque for the first m < 60 modes are shown in Figure [3] for various 
combinations of disk scale height H and viscosity a. As a planet migrates slowly inwards 
from large r, the net Lindblad torque, Ttotab is negative, driving the planet inwards. When 
the planet encounters a disk resonance strong enough such that Ttotai > the planet can 
become trapped. Further migration inwards towards the resonance causes the net Lindblad 
torque to become positive, driving the planet outwards, while outwards migration away from 
the resonance leads to a net negative Lindblad torque, driving the planet inwards, effectively 
trapping the planet where Ttotai — 0. 

In Figure [3l for low viscosity {a = 10~^), we see many locations of effective planet 
trapping for disks with both H/r = 0.1 and 0.05. However for higher viscosity (a = 10^'^) 
we only see trapping for H/r = 0.1, consistent with our estimate (ISTl) . 

Including the effect of the corotation torque Tcr requires a location where ^ T = 
TcR + TiLR + ToLR = to trap the planet. The value of the unsaturated corotation torque in 
an isothermal disk is proportional to the vortensity gradient at the corotation point (GT79), 

dln[K'^/{2J:Q)] 



Tcr oc - 

dmr 



(52) 

''CR 
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and can be larger than the non-resonant net Lindblad torque near sharp density transitions. 
Baruteau fc Massetl (120081 ) also showed that the corotation can contribute a significant torque 
in the presence of entropy gradients in non-barotropic disks. Depending on the sign of the 
density and entropy gradients in the disk at the planet location the corotation torque can 
make it easier or more difficult to trap a planet using the trapped wave disk resonance 
discussed above. 



4. Nonlinear Resonance 



Goodman fc Rafikovl (1200 ll ) and iRafikovl (120021 ) showed that the spiral density waves 



generated by a perturbing planet typically steepen as they propagate and can quickly act to 
shock and damp out the wave. For waves resonant in the disk this would lead to a non-linear 
resonance, in which the primary means of dissipation involves wave shocks. Detailed exami- 
nation of the non-linear response of the disk is left for later work, and may perhaps be best 
approached numerically. However, we note that typically for forced non-linear resonators, 
such as shock tubes, the resonator quality factor, Q, effectively o ur resonant gain, scales 
as the inverse square root of the forced amplitude (JRudenkd l2009l ) . This implies that the 
non-linear resonant gain decreases with with an increase in the planet mass, and that there 
is a cutoff mass above which the resonance can no longer halt Type I migration. 



Discussion 



We have shown that for small planets the linear response of a disk resonance caused by 
an inner reflecting edge can act to trap planets at a few times the inner disk radius. The 
exact location of this trapping depends on the details of the disk parameters and dissipation 
mechanism. For larger planets, the disk will likely undergo a non-linear resonance, which 
typically has non-linear resonant gain decreasing with increasing planet mass. This implies 
that above a certain cutoff mass, the resonant gain due to the non-linear resonance will 
be insufficient to halt migration, and the planets will continue to migrate inwards through 
the inner disk edge. There they will be halted at rhait when the lowest order outer Lindblad 
resonance is located within the inner edge of the disk. The corotation torque is also important 
near these edges, but is susceptible to saturation, unless the viscosity is sufficiently high. If 
unsaturated the corotation torque is mainly active when the planet is very close to the 
disk edge and the density gradient (and thus vortensity gradient) are strongest. However, 
further from the disk edge where the wave trapping resonance is active, both effects may be 
important in determining the net torque. 
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This naturally generates two populations of planets, with halting radius a few times 
the inner disk edge for smaller planets, and with larger planets halted inside the inner disk 
ed ge. This process may help to explain the observations of planets close to solar type stars 



by [Howard et al.l (120 ll[ ) which showed that larger planets were halted with typical periods 
of ~ 2 days, corresponding to rhait ~ 0.03 AU, while smaller planets were halted with typical 
periods of ~ 7 days, rhait ~ 0.07 AU, consistent with inner disk edges with periods of ~ 4 
days corresponding to r^^ ~ 0.05 AU. 

This process may also be important for planets located near other reflec ting edges within 



the disk, such as the edges of magnetic dead zone s fsee e.g. iGammie 119961) . or the edge of a 



gap in the disk opened by a giant planet (see e.g. iLin fc Papaloizou 



197 



i. 



The effect of turbulent density fluctuations on this wave resonance may be to ratchet 
the system through successive "traps" , if the variation is on a scale similar to the separation 
between the Lindblad resonances, such that the relative amplitudes of the Lindblad torques 
are affected. Larger resonant peaks will be more difficult to overcome, however a detailed 
study of the effect of turbulence is beyond the scope of the current work. 
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Fig. 3. — Net Lindblad torque, Ttotai, for the m < 60 modes as a function of planet orbital 
radius Vp. The torque is negative except at the resonant peaks where the gain due to wave 
trapping is high enough for the inner Lindblad torque to overcome the outer Lindblad torque. 
The net Lindblad torque is normalized by Tin = {Mp/M^)'^'Er'^Q'^{H/r)~^ evaluated at the 
inner edge of the disk r = rin. The red lines indicate the differential Lindblad torque without 
considering disk resonance, while the black shows the differential Lindblad torque including 
the effect of disk resonance. When corotation torque is negligible planets can become trapped 
where the net Lindblad torque crosses zero with negative slope. 



